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Abstract 

We evaluate by direct calculation the Lopatinski determinant for ZND detonations 
in Majda's model for reacting flow, and show that on the nonstable (nonnegative real 
part) complex half-plane it has a single zero at the origin of multiplicity one, implying 
stability. Together with results of Zumbrun on the inviscid limit, this recovers the 
result of RoqueJoffre-Vila that viscous detonations of Majda's model also are stable for 
sufficiently small viscosity, for any fixed detonation strength, heat release, and rate of 
reaction. 

1 Introduction 

In this note, we verify by explicit computation the spectral stability in the sense of Er- 
penbeck [Er] of strong ZND detonation wave solutions of Majda's model for reactive gas 
dynamics [M] with a step-type ignition function. 
Consider the inviscid Majda's model 

(1.1) {u + qz)t+ (j) x = °' 

z t + k(p(u)z = 

an analog of the Zeldovich-Von Neumann-Doering (ZND) equations for reactive gas dynam- 
ics, where z, tp £ K, u > 0, q > and k > 0. Here, u is a lumped variable modeling the 
gas-dynamical quantities of density, momentum, and energy, z is mass fraction of reactant, 
q > a coefficient of heat release of the reaction, k > reaction rate, and <p(u) is a simple 
step-type "ignition function" that is assumed to be zero below a certain value ui > and 
one above. 

A strong detonation wave of (jl.ip is a traveling-wave solution 

(1.2) (u,z)(x,t) = (u,z)(x - st), lira. (u,z)(g) = (u±,z±) 
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in the weak, or distributional, sense, smooth except at a single shock discontinuity at 
(without loss of generality) x = 0, known as a "Neumann shock", where u jumps from 
ii* := u(0~) to n(0 + ) as x crosses zero from left to right, and satisfying 

(1.3) z_ = 0, z + = 1, u~ > Ui > u + . 
and 

(1.4) u_ > s > u + > 0. 

Here, analogy of the model for gas dynamics is only for u + > 0, the "physical" range defined 
by Majda. Computing the Rankine-Hugoniot conditions at the shock at x = 0, we find 
that 

, , n , , . u(0+) 2 n(0-) 2 

(1.5) ^(0 + )-^0-)) = ^-^, 

2(0") = ^(0 + ), 

which satisfies 

ti a\ 2 r, 2 o U + U * 

(l.b) u + — lsu + = — zsu* or s = , 

letting U-, q > 0, k > vary. From (ll.4p and the assumption that (u,z) converges as 
x — > ±co, we find further by consideration of the traveling-wave ODE (see Section [2]) that 

(1.7) (u,z)(x) = (u+,z+) for x > 



and also U- = s + J s 2 — 2qs + u\ — 2su + , so that 

(1.8) s < u- < u*, < q < — — , 

2s 

with reaction rate varying in the infinite range < k < +oo. 

That is, we have the standard picture of a strong detonation wave as a shock advancing 
to the right into a quiescent (i.e., nonreacting) constant state with reactant mass fraction 
z = 1, raising u above ignition level Ui, followed by a smooth "reaction tail" in which 
combustion (reaction) occurs, in which z decays exponentially to value z_ = and u to 
s < u_ < -u* as x -> -oo. See [DyZTj [LyZ2l [JLWl IZT1 [Z2] for further details. 

As shown by Erpenbeck [Er], spectral stability of such waves, defined as nonexistence 
of normal modes e xt w(x) with > and A ^ 0, may be determined by examination of a 
certain Lopatinski determinant Dznd{^) whose zeros on $RA > correspond to eigenvalues 
A. Precisely, spectral stability corresponds to the condition 

(D) Dznd(-) has a single zero on {!RA > 0} occurring at A = 0. 
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The Lopatinski condition (D) has been much studied numerically in both one and multi 
dimensions, for the full equations of reactive gas dynamics; see, for example, [KS] and ref- 
erences therein. However, so far as we know, it has up to now not been verified analytically 
for any case. Here, we show by direction calculation that (D) holds for detonations of the 
inviscid Majda model with step-type ignition function. 

Since linear or spectral stability concerns only behavior near values of the profile (u, z), it 
is clear that this result extends also to the slightly more general class of (possibly smooth) 
ignition functions identified by Roquejoffre and Vila |RV] that 4>{u) = for u < U{ and 
4>{u) = 1 for u > u'\ provided that u- > u % > ui > u + . That is, we require in our analysis 
only that 

(1.9) 4>(u) = near u + and 4>(u) = 1 near 

A very interesting open problem is whether a corresponding stability result holds for more 
general ignition functions satisfying only 

(1.10) <M^) = for u < Ui and (j)( u ) > f° r u > u i- 

As discussed in Section [6l these results have implications also for spectral and nonlinear 
stability of viscous detonation waves in the ZND limit as viscosity goes to zero. 

2 Profile solutions 

Observing that 

dtu(x — st) = —su, d x u{x — st) = u d{z{x — st) = —sz, d x z(x — st) = ~z , 
we obtain for a detonation solution (jl.2p - (jl.5p the profile equations 



valid on the regions of smoothness x > and x < 0. Linearizing (|2.1|) about the equilibrium 
solution (n, z) = (n + , 1), we readily find from s > u + that this is a repelling equilibrium, 
and so the only smooth solution of (|2.ip converging to this equilibrium is the constant 
solution (pTT7|) . 

Integrating (|2.ip from — oo to +oo, we obtain 



(2.1) 




sz' + kip(u)z = 0, 



s(u + + qz + ) 




s(u^ + qz-) 



2 



2 



Substituting z + = 1 and z. 



0, we find that 



(2.2) 



2s(u + + q) — u+ = 2su- — u 2 _ 
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Now, integrating (|2.ip from — oo to £, we obtain 



u 2 u 2 



-s{u + qz) + s(u_ + qzJ) = — - + — , 



from which, by (jl.3p and (|2,2p . we obtain 



(2.3) = s ± y^s 2 - 2<?s(l - 2(f)) + «+ ~ 2 su+. 

Recalling is the left-limit of u(£), u* > s and z(Q) = 1, we choose the sign +, yielding in 



particular u- = s + J s 2 — 2qs + u\ — 2su + . Since ip(u) = 1 for £ < in (I2.1j) . z' = (k/s)z, 
hence for £ < 0, 



(2.4) z(0 = Z(0)e = e 

and 



(2.5) u(f ) = s + \]s 2 - 2qs(l - e^) + u\ - 2su + . 

3 Eigenvalue equations 

Let (u,z) be a solution of (II. ip different from (u, z). Subtracting (u, z) from (u,z), we 
obtain for the perturbation variable 

(u,z) := (u,z) - (u,z) 

the perturbation equations 

N U 2 \ 

(« + **)*+ (y-y^ = o 

z t + k(ip(u)z — (p(u)z) = 0. 

Taylor expanding and dropping 0(\(u, z)\ 2 ) terms, we obtain the linearized equations 

ut — su^ + (uu)^ = qk(d(p(u)zu + (p(u)z) 
zt — sz^ = —k(d(p(u)zu + (p(u)z), 

and, finally, the linearized eigenvalue equations 

An — sv! + (uu)' = qk(dtp(u)~zu + ip(u)z) 
Xz — sz = —k(d(p(u)zu + cp(u)z), 

where / denotes d^. 
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Expanding and rearranging, we have 

(II — s)u + (u — s)'ii = (qkd(p(u)z — X)u + qkip(u)z 
—sz' = —kdip(u)zu + (—kifiu) — X)z 

which gives the matrix equation 

u — s \ / u \ 1 ' / qkd(p(u)z — X qkip(u) \ ( u 

—s J \ z J $ \ —kd(p(u)z —kip(u) — X J \ z 

U qkdip(u)z qkip{u) \ _ ^\ ( u 
y \ —kd<p(u)z —kip(u) J 

Setting 



(3.1) A=[ U ~ S ° ),W=( U ),E=( ^/K- ),Z = AW, 



— s J ' y z y ' y —kd(p(u)z —kip(u 

we may write the eigenvalue equations as the first-order ODE system 
(3.2) Z' = (E - XI)W = (E — Xl)A~ x Z = GZ, 

where G = (E-XI)A- 1 . 

4 The Lopatinski determinant 

Following [JLW^ IZ2] , define on $RA > the Lopatinski determinant 
(4.1) D ZND (X) := det(Z-(A,0), X\W] - [AW']), 

where [h] := h(0 + ) — h(0~) and Z~(X,£,) is a bounded exponentially decaying solution of 
(13. 2D . analytic in A and tangent as £ — > — oo to the subspace of exponentially decaying 
solutions of the limiting, constant-coefficient equations Z' = G-Z. By standard asymptotic 
ODE theory (the "gap lemma" [GZJ), Z~ is uniquely determined up to a nonvanishing 
analytic factor not affecting stability. 

Lemma 4.1. The ZND Lopatinski determinant is given (up to nonvanishing analytic factor) 
by 



(4.2) Dznd = ((u* - u+)A + («* - u + - q - qk^)k) ( 



A 



k + X 



„ ik+X) 



(4.3) 



p(0 = 



s 2 - 2gs(l - e ( fc /^) + u\ - 2su A 
X 



s 2 - 2qs(l - e( k / s K) + v? + - 2su. 
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Proof. Combining (|2.4I) . (|2.5p and (J3JJ), we obtain 



A(W] - [A(W )] = A 
= A 



+ (Ar)(o-) 



where we are using W 



G=(E- \I)A~ 



u + — u* \ f u* — s A / (qk)/(u* — s) 







+ 



-s 



Au+ — An* + q'/c 
-jfc 



s 2 - 2gs(l - e( fc /^) + u 2 - 2su H 



fc/a 



k \T 
= ) _ e (*A0£) _ Computing 



—A g/c A / u — s 
-it -A J I -s 



-i 



/ ~ A ~1 k 

u — s s 

*±* 



and setting Z = ( ^ ) , we have 



or 



z' 



\ 



„ - .s .s i ( z i(\,0 
k + X I Z 2 (A,0 



A 



Zj Z^ 



M — S S 

/c + A 



Solving, we obtain 
Z[ 



A ry Qk ( fc + A ) f 

Z\ e « S 

s 2 - 2gs(l - e ( fc / s )5) + u 2 , - 2su+ s 



Z9 = e s 6 . 



Setting now P(£) 



A 



s 2 - 2qs{\ - e( fc / s )5) + u\ - 2su+ 



and := e s > we 
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have 
(4.4) 



k\ s 2 — 2qs + u 2 , — 2su+ 



x In 



s 2 - 2qs{\ - e( fc / s )?) + u 2 - 2su + + J s 2 -2qs + u\-2 



V / 2qse ( - k / s ^ 



which gives 



e-Z-oo-PCO* 



< 1 for 5RA > 0, and hence, integrating by parts, 
Zi(A,0) = e -f-~ P{s)ds ( f e^ P{s)ds Q(y)dy 

\J—oo ' 

ti P ^ ds Q(y)dy 



(4.5) 



o 

e 



fe + A 



-oo J — oo 



:i-a^, 



A; + A 

where * is as in (J33J). By (|3~Tj) . (Oj) . and fl33]), 

DzndW = det(Z-(A,0), A[W] +^(W')(0-)) 

_ , / ^i(A, 0) Au + - Au* + g/c 
- det V^ 2 (A,0) -fc 

^4 g\ = — fcZi + A-u* — \u + — qk 

_ qk 2 (l - A*) + (Au* - Au+ - qk)(k + A) 
~ k + A 

A 



((it* - u + )A + (it* - u+ — q — qk^)k) 



k + X 



□ 



5 Verification of spectral stability 

Theorem 5.1. For a step-type ignition function, or, more generally, any ignition function 
satisfying flL9j) , Dznd has a single zero of multiplicity one on {3?A > 0}, located at A = 0; 
that is, the Lopatiski condition (D) is satisfied for all ZND detonations of Majda's model, 
independent of the choice of q > 0, k > 0, or u±. 
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Proof. By ([0|), we have for KA > 



o (k+ntx) 
-°° A /s 2 - 2qs(l - e( k / s )y) + ul - 2su + 



fO e (k/s)y 

< / , =<fr 

's 2 - 2gs(l - e( fc /%) + u 2_ _ 2s-u + 



-|- (^/s 2 - 2gs(l - e( fc / s )f ) + u\ - 2sii+) cfy 
fy/ s 2 + u+ — 2su + — \j s 1 — 2qs + u 2 + — 2su^ 



s 2 — 2qs + u\ — 2su + 



giving 



(5.1) 



K(u* — u + — q — qk^) > u* — u + — q — qk\^\ 



> it* — u + — q — s + u + + \ s 2 — 2qs + u + — 2su 



s — q + y s 2 — 2qs + u + — 2su + 



s — u + — q + \/ s 2 — 2qs + u 2 + — 2su^ 
By < u + < s and 2qs < (s — u + ) 2 < s 2 , 



K(u* — u + — q — qk^>) > 2qs — q + \ / s 2 — 2qs + u 2 + — 2su 



( 5 - 2 ) > - q + yjs 2 - 2qs + u\ - 2su + 

> 0. 

In particular, 

(5.3) u* - u + - q - qk$> + for KA > 0. 

Combining fj4.6|) and ()5.3|) . we obtain the result. □ 

6 Viscous stability and a result of Roquejoffre— Vila 

Our results have implications also for stability of "viscous" detonation waves, i.e., smooth 
analogs of traveling waves (jl.2p satisfying the "viscous" or parabolic regularization of (jl.ip : 

^g-Q (u + gz)t+^y^ = eqz xx + eu xx , 

z t + kip(u)z = €Z XX , 
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with e > 0. It is shown in [Z2J that stability of viscous detonation waves in the ZND 
limit as e — > with other parameters held fixed is equivalent to Lopatinski stability of 
the limiting ZND detonation (II, 2D together with Evans (equivalently, linearized) stability 
of the viscous Burgers shock corresponding to the Neumann shock contained in the ZND 
detonation profile, as is well-know to hold by properties of scalar traveling-waves (see, e.g., 

[EI]). 

Thus, our results together with those of jZ2] yield the result for (16. lh of spectral stability 
of viscous detonation waves in the ZND limit, similar to an earlier result of Roquejoffre- 
Vila [RV| for the corresponding equations with regularization eu xx in the u equation alone 
and applying to the same class of ignition functions (jl.9p . A very interesting open problem 
would be to extend our results to the more general class of ignition functions fll . lOf) , which 
would give new information for the viscous stability problem as well. 

We remark that results of [LRTZ] (for Majda's model) and |TZ] (for the physical reac- 
tive Navier-Stokes equations) show that spectral stability of viscous detonation waves in 
the (Evans function) sense of [Z2| implies linearized and nonlinear orbital stability, hence 
stability of strong viscous detonation waves reduces for viscosity sufficiently small to a study 
of spectral ZND stability as carried out here. 
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